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Abstract. - We theoretically study sedimentation-diffusion equilibrium of dilute binary, 
ternary, and polydisperse mixtures of colloidal particles with different buoyant masses and/or 
charges. We focus on the low-salt regime, where the entropy of the screening ions drives 
spontaneous charge separation and the formation of an inhomogeneous macroscopic electric 
field. The resulting electric force lifts the colloids against gravity, yielding highly nonbarometric 
and even nonmonotonic colloidal density profiles. The most profound effect is the phenomenon 
of segregation into layers of colloids with equal mass-per-charge, including the possibility that 
heavy colloidal species float onto lighter ones. 



Introduction. - Sedimentation in a suspension of colloidal particles is of profound funda- 
mental importance, and has been studied in detail for a long time. For instance, in 1910 Perrin 
determined the Boltzmann constant fc^ (and from this Avogadro's number) by comparing the 
measured equilibrium density profile of a dilute suspension with the theoretically predicted 
barometric law [1], and more recently the full hard-sphere equation of state was determined 
accurately from a single measurement of the sedimentation profile of a dense suspension of 
colloidal hard spheres [2,3]. In rather dilute suspensions of charged colloids, however, strong 
deviations from the barometric distribution have recently been theoretically predicted [4-9], 
experimentally observed [10-13], and simulated [14], at least in the regime of extremely low 
salinity. The most striking phenomenon is that the distribution of colloids extends to much 
higher altitudes than is to be expected on the basis of their buoyant mass. The force that 
lifts the colloids against gravity is provided by an electric field that is induced by sponta- 
neous charge separation over macroscopic distances. The corresponding cost of electrostatic 
and gravitational energy is more than compensated by the gain of translational entropy of 
the salt ions, which are distributed much more homogeneously than when they would have 
had to follow the colloids to the bottom part of the suspension. The same entropy-induced 
lifting force was recently found to be responsible for a so-called Brazil-nut effect in binary mix- 
tures of charged colloidal particles under low-salt conditions: the heavy particles can reside 
at higher altitudes than the lighter ones [15]. In this Letter we extend the theoretical study 
of sedimentation of charged colloids to binary, ternary, and polydisperse mixtures, for which 
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we calculate colloidal density profiles and the entropy-induced electric field selfconsistently 
within a Poisson-Boltzniann like density functional theory. 

Theory. - We consider an n-component suspension of negatively charged colloidal spheres 
in a salt solution in the Earth's gravity field. We label the colloidal species by i = 1, 2, • • ■ , n, 
and denote the electric charge, diameter, and buoyant mass for species i by —ZiC, <Ti, and m^, 
respectively, where e is the proton charge. We imagine the suspension to be in thermal and 
osmotic contact with a reservoir that contains monovalent ions at a concentration 2ps and a 
solvent with dielectric constant e at temperature T . The ions are assumed to be massless point 
particles with charge ±e, and the solvent mass density is taken into account by considering the 
buoyant instead of the actual colloidal masses according to Archimedes. The volume of the 
suspension is F = AH , with A the (macroscopic) horizontal area and H the vertical height of 
the solvent meniscus above the bottom of the system. The gravitational acceleration g is in 
the negative vertical direction. We are interested in the equilibrium colloidal density profiles 
Pi{x) as a function of the altitude x above the bottom at x = and below the meniscus at 
X — H. The ion density profiles are denoted by p±{x). In order to calculate these profiles we 
employ the framework of density functional theory [16], where the equilibrium density profiles 
follow from the minimisation of a grand potential functional U,[{pi\, p±\ with respect to all 
the profiles. Here we employ the mean-field free energy functional that consists of entropic, 
gravitational, and electrostatic contributions, 
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where we introduced the fugacity (or activity) and the gravitational length Li = ksT /rriig 
of species «, the Bjerrum length Xb — e^/^ksT, and the local charge density Q{x) — p+{x) — 
P-ix) — '^"^i ZiPi(x). The electrostatic term follows from an in-plane integration of the three- 
dimensional Coulomb law, 27r dRR/ y^\x — a:'p -f i?^ = — 27r|a; — x'\ with R the in-plane 
distance, where we ignored an irrelevant integration constant. The functional only couples 
the total charge density at different heights in a mean-field fashion, i.e. electric double layers 
and many other correlation effect are not taken into account. In fact one easily checks that the 
thermodynamics of the system reduces in the absence of gravity to an (n-|-2)-component ideal- 
gas mixture, since then Q{x) = because of translational invariance and charge neutrality. In 
the presence of gravity, however, interesting structures already appear at this relatively low 
level of sophistication. The Euler-Lagrange equations 5Q./5pa{x) = and Sn/Spi{x) = that 
must be satisfied by the the equilibrium density profiles can be cast in the form 

P±{x) = Ps exp[T(t)ix)]; (2) 
Pi{x) = fli exp[-a;/Li + Zi(/)(a;)]; (i = 1, 2, • • • , n) (3) 

where (j){x) = —2t:Xb dx'Q{x')\x — x'\ is the dimensionless electrostatic potential gauged 
such that it is zero in the reservoir where Q{x) = 0. It turns out to be convenient to rewrite 
it in differential form as the Poisson-Boltzmann equation 

= -47rAi3Q(x) = K^sinh0(x) -|-47rAB^Z,p,(x). (4) 

i—1 
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Here we used Eq.Q, and we defined the reservoir screening constant — SttXbPs- Together 
with the boundary conditions 0'(O) = 4''iH) — (where a prime denotes a derivative w.r.t. x) 
the Eqs.lO and Q form a cfosed set of n + 1 equations that can in principle be solved to yield 
Pi{x) and 0(2;) for a given thermodynamic state determined by the reservoir characteristics 
K and Ab, and by the colloidal fugacities a^, gravitational lengths Li, and charges Zi. In an 
experimental situation, however, the total packing fraction f]i = {tt /Q)Niaf /V of species i 
is usually fixed instead of the fugacity a^, with Ni the number of colloids of species i. This 
conversion can easily be accomplished here by regarding the fugacities as normalization 
constants that take values such that fji = (l/H) J^^ dxrii{x), with rii{x) = {Tr/6)pi{x)af the 
local packing fraction of species i. Note that the hard-core diameters ai do not appear at 
all in the Euler-Lagrange equations, since the functional Q ignores the hard-core part of the 
(direct) correlations; we use for all colloidal species that (Ti = a = 150nm only to be able to 
convert densities to physically reasonable packing fractions. Once the profiles (f>{x) and Pi{x) 
are known, then p±{x) follows from Eq.Q, and from this the total charge density Q{x). The 
magnitude of the electric field is given by E{x) — ksTcj)' [x) / e, where the prime denotes a 
derivative w.r.t. the height x. 

Unfortunately one cannot solve Eqs.© and Q analytically. However, its one-dimensional 
character allows for a rather straightforward numerical solution on a grid {xfc} of heights, 
although some care must be taken in dealing with the widely different length scales H (say of 
the order of centimeters) and (at most of the order of microns) in realistic cases. In all 
our calculations we use a non-equidistant grid that we adapt to the situation at hand, making 
sure that it is fine enough to exclude numerical artifacts. Typically we use 600 grid points for 
a system with a meniscus ai H = 20cm. Our iterative scheme to solve Eqs.Q and Q on the 
grid takes typically a few seconds on a desktop PC for n — 1, and a minute for n = 21. 

Before discussing our numerical results, we wish to point out that the present theory 
reduces, for n = 1, to the treatment of Ref. [8]. For instance, the three regimes for the colloidal 
density profile (barometric, linear, and exponential with a large decay length) follow within 
the assumption of local charge neutrality Q{x) = 0, i.e. —s\nh.if){x) = Zipi{x)/2ps = yi(x). 
In the linear regime, where 1/Zi < yi{x) < 1, the potential is then from Eq.l^J linear in x, 
with a slope (f>'{x) = 1/ZiLi that is proportional to the electric field that lifts the colloids 
against gravity [8] . This one-component result already hints at an important complication for 
mixtures (n > 2): since generally ZtLi ^ ZjLj a single electric field strength cannot lift all the 
colloidal species simultaneously. On this basis one could therefore already expect segregation 
in mixtures, such that colloids with the same value ZiLi (i.e. the same charge-per-mass) are 
found at the same height. This is indeed what our numerical results will show below. 

Segregation in binary and ternary mixtures. - We start our numerical investigation with 
a class of binary mixtures {n = 2) of equal-sized light {i = 1) and heavy (i = 2) colloidal 
particles of various charge ratios. We choose the system parameters identical to those of the 
Monte Carlo simulations of figure 3 of Ref. [15]: a ~ 150nm, H — lOOOcr, Xb = cr/128, 
Li = 1^2 = lOcr, ^1 = 772 = (7r/6) X 10^"*, Zi = 15, and Z2 varies between 15 and 45. The 
simulated system does not contain added salt but only counterions, which we can represent 
within our theory by the extremely low reservoir salt concentration ps — InM, which we 
checked to be low enough to be in the zero-added-salt limit as regards the colloidal profiles. 
In Fig. 1(a) we show, for Z2 = 45, the density profiles of the two colloidal species, as well as 
that of the counterions (and the colons) in the inset. We observe profound colloidal segregation 
into two layers, with the heavy species floating on top of the lighter ones. The counterions 
are seen to be distributed throughout the whole volume, i.e. much more homogeneously than 
when all colloids would have been barometrically distributed in a thin layer of thickness Li 
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just above the bottom (since in that case the net ion charge would have been located in that 
same thin layer). The resulting gain of ion entropy is the driving force for the formation 
of the electric field that pushes the heavy (highly charged) colloids to high altitudes against 
gravity [8,14]. In Fig. 1(b) we show the mean height hi of species i, defined as 



_ dx X pi{x) 



hi = 



jf dx pi{x) 



(5) 



as a function of Z^. We replot the Monte Carlo simulation results of figure 3 of Ref. [15] 
(symbols) , together with the predictions for hi that follow from the present theory (continuous 
curves). Given that there is not a single fit parameter involved, the agreement is remarkable, 
certainly when compared with the theoretical analysis on the basis exponentially decaying 
density profiles as in Ref. [15]. Fig.l shows that the heavy particles are on top of the lighter 
ones, /i2 > hi, provided ZijZx ~ 1.6, a phenomenon that was termed the "colloidal Brazil 
nut effect" in Ref. [15]. For barometric profiles one would find that hi = Li, but we see in 
all cases that hi » Li due to the lift effect of the induced electric field. The good agreement 
between theory and simulation in Fig. 1(b) also indicates that hard-core effects (which arc 
taken into account in the simulations but not in the theory) are not so relevant, at least not 
in the parameter regime studied here. 
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Figure 1 - (a) Colloidal density profiles and counterion distribution (inset) for a deionised binary 
mixture of equal-sized light (i = 1) and heavy {i = 2) colloidal spheres (diameter a = 150nm), 
with meniscus height H = lOOOcr, gravitational lengths ii = lOcr and L2 ~ 20cr/3, average packing 
fractions 771 = ri2 = (7i"/6) x 10^*, colloidal charges Zi = 15 and Z2 = 45, and Bjcrrum length 
Xb = cr/128. (b) Mean height hi of the two species for the same system as in (a) except now for a 
range of Z2, showing a density inversion (/i2 > hi, the heavy particles floating on top of light ones) 
for Z2/Z1 > 1.6. The curves are predictions of the present theory, the symbols are simulation data 
of Ref. [15]. 



The next system we study is a ternary system (n = 3) in a solvent characterised by 
Xb = 2.3nm (ethanol at room temperature) with a reservoir salt concentration ps = lOfiM 
and meniscus height H = 20cm. The colloidal charges and gravitational lengths are Zi = 
(1000, 250, 125) and Li = (1, 2, l)mm, respectively, and the system is equimolar with fji = 
0.005 for all three species (recall that all species have the same diameter <t = 150nm). Fig. 2 
shows the density profiles as predicted by the present theory in (a), as well as the electric field 
in (b). We find almost perfect segregation into three layers, such that the mean heights satisfy 
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< h2 < hi. Note that this ordering coincides with the ordering < Z2L2 < ZiLi, and 

not with the corresponding ordering of Li which one would expect on the basis of a barometric 
distribution. In other words, this system segregates according to mass-per-charge instead of 
the more usual ordering according to mass: the colloids with the largest mass-per-charge are 
found at the bottom. In the present case this implies that the lightest colloids (species 2) 
are found in a layer in between the equally heavy species 1 and 3. On the basis of the one- 
component theory of Ref. [8] one would expect a linearly decaying density profile of species i in 
the layer where species i is dominant, as well as a linear electrostatic potential. Denoting the 
derivative w.r.t. a; by a prime, the one-component theory predicts = irpsU^ /iZfLi and 

4>' {x) = l/ZiLi in the layer with species i, corresponding to an electric field strength rriig/Zie 
in this layer. These values for the density gradients and the electric field are indicated in 
Fig. 2(a) and (b), respectively, and are in good agreement with the numerical results. The 
result for the electric field can also be easily obtained analytically for a mixture provided one 
assumes that segregation takes place, such that pi [x] takes a maximum at some height x* in the 
layer of species i. From p'i{x*) = one obtains from Eq.® that (t>'{x*) = 1/ZiLi, in agreement 
with the numerical results. The inset of Fig. 2(b) shows the ratio of the total charge density 
Q{x) and the ion charge density 2psSmh(j){x), which is such that |Q(a;)/2/3<; sinh(/)(a::)| ^ 1 
for all X except close to a; = and x — H where it is ~ 1. This indicates that the system, 
with its inhomogeneous electric field, is yet essentially locally charge neutral (but not exactly, 
and not at all at the boundaries), suggesting that a description on the basis of hydrostatic 
equilibrium within a local density approximation should be rather accurate [17]. 
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Figure 2 - Density profiles (a) and the electric field profile (b) of an equimolar ternary mixture of 
equally-sized colloids with charges Zi = (1000,250,125) and gravitational lengths Li = (1,2, l)mm 
with total packing fractions fji — 0.005 in a 20cm ethanol suspension at room temperature, with 
monovalent ionic strength ps — lOpM in the reservoir. Almost complete segregation takes place into 
well-defined layers of pure components, with a slope r]i{x) — ivpsa^ /SZ^ Li and a constant electric field 
rriig/eZi in the layer with component i, as denoted by the dashed lines in (a) and (b), respectively. 
The inset in (b) shows a dimensionless measure for local charge neutrality (see main text). 



Polydisperse Mixtures. - We now extend our study to polydisperse mixtures, where 
one could expect segregation into many layers on the basis of the results for two and three 
components. We mimic the polydispersity by considering a system of n = 21 components, 
with Zi distributed as a Gaussian with average of 250 and a standard deviation of 50. This 
distribution is shown in the inset of Fig. 3(a), where the vertical axis {fji) is proportional to 
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the relative frequency of species i in the sample. We consider two distributions for Lf. (A) 
Li — 2mm for all species, and (B) Li — 2 y. (250/Zi)'^/^mm. Case B mimics the situation 
for spheres of different size but the same mass density and surface charge density, such that 
Zi is proportional to the surface area and Li to the inverse volume of species i, i.e. LfZf 
is a constant independent if i. The density profiles, numerically obtained by solving the 
Eqs.© and Q for iJ = 20cm, Xb = 2.3nm, ps = 3pM, and fjtot = T.'Him = 0.005, are 
shown in Fig. 3 for all 21 components. Case A shows profound lifting and layering, where the 
ordering is again determined by mass-per-charge as illustrated by the three dashed curves (for 
Zi = 250, 279, and 299) showing that the colloids in the high-charge wing of the distribution 
reside at high altitudes. Fig. 3(b) shows the density profiles for case B, which does exhibit 
lifting, but hardly any layering, and no density inversion at all. This is completely consistent 
with the picture that the ordering is determined by ZiLi, which for case B is such that 
ZiLi/ZjLj = ^ Zj I Zi , i.e. the highly-charged particle are expected at the bottom while the 
relative spread in ZiLi is relatively small compared to case A, where ZiLi/ZjLj = Zi/Zj. 
The inset of (b) shows the total packing fraction profiles ritot{x) = Yl'iLiVii^) of both case 
A and B together with the one-component profile {n — 1) with Zi = 250 and Li — 2mm 
at fji — 0.005. Perhaps surprisingly there is hardly any distinction between the pure system 
and case B, whereas there is a small difference with case A. These profiles show that the 
main distinction between these polydisperse systems and the underlying one-component one 
concerns the layering phenomenon (provided ZiLi varies sufficiently for all the species), and 
not the total distribution of the colloids. Perhaps this fractionation effect could be exploited 
experimentally to purify a polydisperse mixture. 
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Figure 3 - Density profiles of a 21-component colloidal suspension with a Gaussian charge distribution 
as illustrated in the inset of (a), at total packing fraction 0.005, reservoir salt concentration ps — 3/iM, 
and Bjerrum length Xb ~ 2.3nm. The gravitational lengths are as in case A (see main text) in (a), 
and case B in (b). In (b) the solid curves are for Zi > 250, and the dashed ones for Zi < 250. 
The inset of (b) shows the total packing fraction profile for case A and B, together with that of the 
underlying one-component system (see main text). 

Conclusions and discussion. - We have studied sedimentation equilibrium of n-compo- 
nent systems of charged colloidal particles at low salinity within a Poisson-Boltzmann-like 
density functional. The Euler-Lagrange equations that describe the minimisation of the func- 
tional are solved numerically on a one dimensional grid of heights for a number of system 
parameters. For n = 1 the theory reduces to the one-component studies as presented in 
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Ref. [8,14], and for n = 2 we quantitatively reproduce the simulation results of Ref. [15], 
where density inversion was found such that the heavier colloids float on top of the lighter 
ones. This effect is caused by a self-consistent electric field that lifts the higher charged (heavy) 
particles to higher altitudes than the lower charged (lighter) colloids. We show that the layer- 
ing of the colloids according to mass-pcr- charge can persist for ternary (n = 3) as well as for 
poly disperse (here n = 21) mixtures. Given the good account that the present theory gives for 
simulations [14] and experiments [11-13] of one-component systems, and for the simulations of 
binary systems [15], it is tempting to argue that the theory is also (qualitatively) reliable for 
ternary or polydisperse systems, i.e. the predicted segregation and layering should be experi- 
mentally observable. One should bare in mind, however, that the present theory ignores the 
hard-core of the colloidal particles, and is therefore expected to break down at higher packing 
fractions. It also ignores effects due to charge renormalisation, which becomes relevant when 
ZiXsjoi » 1. Work on extending the theory in these directions is in progress [17, 18]. 

All results presented in this paper were obtained with the zero-field boundary conditions 
(/)'(0) = (p'{H) = 0. We checked explicitly, however, that other boundary conditions that 
respect global charge neutraUty, such as ^'(0) = (j)'{H) — eEextlksT (describing a suspension 
in a homogeneous external electric field Eext) or (f){Q) = <j){H) and = 4>'{H) (describing 
a short-circuited bottom and meniscus) give indistinguishable density profiles, except in two 
layers of thickness ~ <C lO/itm in the vicinity of the bottom and the meniscus. This 
insensitivity to the boundary conditions is not surprising in the light of the fact that the 
whole phenomenology in these systems is driven by the entropy of the microscopic ions, i.e. 
a bulk contribution to the grand potential that should dominate any boundary (surface) 
contribution. 
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